6, 2008 7:41 WSPC/INSTRUCTION FILE K'Vl 



International Journal of Modern Physics E 
© World Scientific Publishing Company 



TERMINATING STATES AS A UNIQUE LABORATORY FOR 
TESTING NUCLEAR ENERGY DENSITY FUNCTIONAL 

M. ZALEWSKI* 

Institute of Theoretical Physics, University of Warsaw, 
ul. Hoza 69, 00-681 Warsaw, Poland 

W. SATULAt 

Institute of Theoretical Physics, University of Warsaw, 
ul. Hoza 69, 00-681 Warsaw, Poland, 

Received (received date ) 
Revised (revised date ) 

Systematic calculations of favored signature maximum-spin I ma x and unfavored signa- 
ture Imax — 1 terminating states for [/^y 2 ] and [djj/j/jyij' 1 '] configurations (n denotes 
number of valence particles) in A~44 mass region are presented. Following the result 
of Zduhczuk et al., Phys. Rev. C71 (2005) 024305 the calculations are performed us- 
ing Skyrme energy density functional with empirical Landau parameters and slightly 
reduced spin-orbit strength. The aim is to identify and phenomenologically restore rota- 
tional symmetry broken by the Skyrme-Hartree-Fock solutions. In particular, it is shown 
that correlation energy due to symmetry restoration is absolutely crucial in order to re- 
produce energy splitting E{I max ) — E(I ma x — 1) in [/™/ 2 ] configurations but is relatively 
less important for [^A/"/^ 1 ] configurations. 

1. Introduction 

Structure of heavy nuclei may be described using two major theoretical approaches: 
the nuclear shell-model^ (SM) and the nuclear density functional theory^ (DFT). 
The main advantage of SM is proper treatment of many-body correlations among 
limited number of valence particles. The resulting wave functions are eigenstates 
to symmetry invariants of the effective Hamiltonian. In the language of DFT, nu- 
cleus is treated as A-hody composite object described by one-body densities and 
currents. According to the Hohenberg-Kohn-Sham (HKS) theorem^, this method 
is formally exact. However, since the HKS theorem does not provide any method 
of constructing such an exact functional standard procedures employed in nuclear 
physics attempt to construct the functional in a systematic way guided by basic 
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symmetry requirements^. A particu lar cas e of the nuclear EDF is the functional in- 
spired by well studied Skyrme model EEE (S_EDF) which, in the isoscalar-isovector 
t = 0, 1 representation, takes the following form: 



g Skyrme = ^ j ^ 
t=0,l 



n ( r\r) + n^°\r) , (1) 



where 

H[ TE) (r) = C p t p\ + C? p Pt A Pt + dptn + C J t J$ + C7 J Pt V ■ J u (2) 

H t (TO) (r) - CIs* + C t As s t As t + Cjs t ■ T t + C\f t + C7 3 s t ■ (V x j t ). (3) 

The functional H. is uniquely expressed as a bilinear form of time-even (te) p, r, J 
and time-odd (to) s,T,j local densities, currents, and by their derivatives. Exact 
expressions linking auxiliary Skyrme parameters Xi, ij, i = 0, 1, 2, 3 and W, a with 
coupling constants C can be found, for example, in Ref. ^. 

The coupling constants of the S-EDF (or Skyrme force (SF)) are adjusted to 
global properties of nuclear matter and to selected experimental observables in order 
to account for basic nuclear structure properties. Since there exist a multitude of 
observables and no real consensus about which of these are to be selected in a 
unique manner, there exist a multitude of different SF parameterizations. One of 
the basic problems in adjusting the force parameters is related to the fact that the 
single-particle (SP) states, that are so crucial for high accuracy calculations, are 
in general coupled to collective motion and therefore difficult to determine. In this 
context, any dataset representing SP motion is an invaluable source of information 
that can be used for a rigorous test and subsequent fine-tuning of the parameters. 

Our long standing experience tells us that terminating states which are 

maximum-spin states with in g iven SP configuration are one of the best examples 

17161 

of unperturbed SP motion LL1S1 . This conclusion was recently confirmed by a com- 
parative study between state of the art SM and Skyrme-Hartree-Fock calculation^ 
performed for maximum spin I rnax states terminating within fy 2 and d^^f^ 1 
configurations in N>Z, A^AA nuclei. Indeed, for these cases the energy difference 

AE = E[d-^+ l ] Imax E[f? /2 ] Imax , (4) 

calculated using Skyrme-Hartree-Fock (SHF) approach with modified spin-fields 
and spin-orbit (SO) strength 351 follows very closely the results of SM calculations. 
This study revealed simultaneously that in N=Z nuclei the extreme SP picture 
breaks down due to spontaneous breaking of isobaric symmetry. 

The aim of this work is to explore the structural simplicity of maximum-spin 
Imax (favored signature) as well as for I ma x — 1 (unfavored signature) terminating 
states in order to further constrain parameters of the S-LEDF as well as in order 
to identify, quantitatively evaluate and subsequently restore in a phenomenologi- 
cal way broken symmetries inherently obscuring the SHF treatment. In particular, 
we will show that the energy difference of Eq. (JU) between I max configurations 
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an d ^5/2 m ^4~16 mass region and between I max configurations /^ffg/s; 1 
and g^ 2 in A^80 mass region follows closely the trend found in A^AA mass region 

in Refs. with respect to the isoscalar effective mass scaled isoscalar strength of 
SO term and that this result nicely correlates with simple Nilsson model predictions 
concerning N=Z=8 and N=Z=A0 magic gaps. In the second part we will demon- 
strate that unfavored signature I ma x — 1 terminating states in A~AA mass region 
manifestly violate rotational invariance. However, unlike in most other cases where 
the spontaneous breaking of rotational symmetry (SSB) leads to the occurrence of 
deformation what allows to incorporate substantial fraction of many-body correla- 
tions into a single deformed Slater determinant, see Ref. ^and refs. quoted therein, 
the effect discussed here occurs at quasi-spherical shape. It appears that there 
is a fundamental difference between these two situations. In deformed nuclei the 
SSB mechanism works constructively . The violated symmetry can be subsequently 
approximately restored within independent particle model based on cranking ap- 
proximation what leads naturally to emergence of collective rotational states. In 
the case of I m ax — 1 terminating states the SSB mechanism works destructively and 
one need to go beyond mean-field and perform configuration mixing calculations to 
restore the symmetry. However, since the number of participating (dominant) con- 
figurations is very limited such configuration mixing can be performed analytically, 
as shown in detail in Sect. [31 



2. Isoscalar spin-orbit strength 

In Refs. it was shown that the set of terminating states in N^Z, A^AA mass 
region provides reliable constraints on time-odd fields and the strength of isoscalar 
SO interaction of the S-EDF. In particular, it was demonstrated that constrain- 
ing coupling constants of time odd spin-fields to the empirical spin-isospin Landau 
parameters leads to unification of predictions for AE of Eq. Q for such different 
Skyrme parameterizations like SLy4 and SLy5 031 skO 01 SII1 d, and SkXc EGA 
with mean deviation between calculated and experimental values Ai?^500keV i.e. 
at the level of ~10%. Let us stress that our high-spin estimate of Landau pa- 
rameters is fairly consistent with the data extracted based on completely di fferent 
experimental input like giant resonances, beta decays, or moments of inertia H^EII 
The remaining mean discrepancy AE = AE exp — AEth between theory and 
experiment can be further reduced by slightly reducing SO strength. Indeed, the 
quantity AE of Eq. {4| is governed by the size of N = Z = 20 magic gap Ae2o 
which, in turn, strongly depends on the S O st rength. It can be nicely demonstrated 
within the spherical Nilsson Hamiltonian 01 where: 

Ae2o = fiLUQ(l — 6k — 2k/i). (5) 

For light nuclei \i ~ i.e. flat bottom and surface effect do not play important role. 
The width of the potential well Hujq determines global energy scale in low energy 
nuclear physics. It is rather well constrained by data and even small variations 
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can spoil in general good agreement between theory and experiment, especially in 
heavy nuclei. Hence, the SO term k plays indeed a dominant role in magnitude of 
the magic gap Ae2o- 
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Fig. 1. Mean value AE as a function of the isoscalar strength Wq (filled diamonds) and Wo 
(open cir cles ') for different parameterizations of the SF. See text for further discussion. Taken 
from Ref . OH. 



In Refs.Q@it was shown, that AE in A^44 nicely correlates with isoscalar SO 
strength. However, due to non-local momentum dependent terms, such a correlation 
cannot be done at the level of bare SO strength Wo, see black dots in Fig. [21 but must 
be performed at the level of the so called asymptotically equivalent representation: 

*V = ]^^' (6) 

In this representation free particles in the infinity r — > oo acquire not effective m* 
but bare mass m while the SO potential takes the form: 

m*(r) ( 1 1 1 - 

V LS (q,r) « W -p' (r) ± W^p'^r) Is. (7) 

m y r r J 

Note, that the true SO strength felt by a nucleon is W * = ^W and not W . The 
correlation between the effective-mass-scaled isoscalar strength is Wq and AE is 
now evident, see diamonds in FigO All forces giving similar AE including SLy4, 
Sly5, SIII, SkO and SkXc have also similar W^ (~ 135 ± 10 MeV fm 5 ). For SkpSS 
MSkl or SkM* EQl w hich give unacceptably large AE, Wq is also considerably 
larger. 

To further investigate correlation between AE and the SO strength Wq we have 
performed similar calculations for terminating states AE — E(Jp^} 2 d^f^]i maal ) — 
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Fig. 2. Left figure shows mean value of AE — E re f in A~44 nuclei in comparison to analogical 
quantity calculated for AE = ^df^^g^^]) — E([g^ 2 ]) configurations in /1~80 nuclei. Right 

figure shows AE — E re f calculated for AE = EQp^Ld^'f^]) ~ ^([^5/2]) conn g ura tions m A~W 
nuclei. See text for further details. 



r-i 0.5 
> 

0.0 



« -0.5 
-1.0 



SLy5 



SLy4 



SIII 
SkO 



*5/2J LS9/2J 

• WW" 1 



Q • 

O 



SkXc 



SkP 



O 



SkM* 

O 



MSkl 

O 



130 140 150 160 170 

Wo [MeV fm 5 ] 



Fig. 3. Comparison of AE — E re f for two possible particle-hole excitations [/, 
\Pl/ 2 9gt^] through the magic gap 40. 
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£([dg/ 2 ] i maa ) in A^IQ mass region and for two different terminating configura- 
tions° AE = E([f- 1 2 g;+ 1 } Ima J - E({g% /2 } Ima J and AE = E(\p-^+ 1 ] Imax ) - 
E{[g§/2[i max ) in A~80 mass region. The results are depicted in Figs. and [3J Since 
there is no data for terminating states in these nuclei we compare purely theoretical 
trends by showing mean value AE — E re f (averaging over different N^Z nuclei) 
where E re f denotes mean value of AE calculated using SkO, SIII and Sly4 forces 
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i.e. those which perform best in A~A4 region. 

Note, that dependence of AE — E re f on Wq for [p^e^ 1 ] and [f^/\9g1^] con- 
figurations is very similar to each other and to A~44 shown in Fig. O For 
\Pi/ 2 9g/2 ] states one observes completely different pattern, see Fig. [3J The reason 
for that becomes clear when one looks at the hierarchy of different contributions to 
the corresponding energy gaps emerging within Nilsson model: 

Afe^^oa-^), (8) 

Ae^ 2 z Z d 2 /2) = k»o(l - 6k - 2 Kf x), (9) 

Aefe^ 5 /2) = tkJoO- -8k- 3Kfi), (10) 

AcnTz'^ = ^o(l - 6k - 13k/*). (11) 

The value of /x vary from zero in light nuclei where the surface and flat bottom 
effects are small to the pseudo-spin limit /i ~ 1/2 in heavy- nuclei. In (|HJ), (© and 
pop contribution to Ae from flat bottom and surface effect is small in comparison 
with contribution from SO interaction. The ratio of these two terms (assuming 
fi = 1/2) for |9]) and (fT0|) cases is equal: 

f«f«0.17, (12) 

and for ([8]) it is equal zero. In case of ph excitation from pi/2 to gg/2 sub-shell this 
ratio is: 

(13) 

and one may anticipate qualitative change in physical scenario what indeed takes 
place as shown in Fig. [3] 



3« Imax 1 states 

Unfavored signature [/"^lu.-i terminating states can be obtained within mean- 
field theory by changing either the signature of valence neutron or the signature 
of valence proton. These two independent Slater determinants will be labeled \v) 
and \ir), respectively. Energy difference, between these states and I m ax solution 
calculated using modified SkO force is shown in the left hand part of Fig. [4] Note 
that the SHF values are in complete disagreement both with experimental data and 
SM results. Similar result holds for SLy4 force, see Ref. ^H. 

The major source of disagreement is related to spontaneous violation of rota- 
tional symmetry by these quasi-spherical SHF solutions. Indeed, simple m-scheme 
counting shows that both I max and I max — 1 representations are single-folded within 
the fy 2 Hilbert space irrespective of number of valence particles n. The SHF solu- 
tions are eigenstates of an angular moment projection with K — I max — 1. Hence, 
they are mixtures of " spurious" \I max ,Imax - 1) and physical \I max - 1, W - 1) 
states: 



\Imax,Imax - 1) = a\ v ) + » 



(14) 
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\Imax — 1, Imax - 1) = -b\v) + a\ic) . (15) 

The symmetry restoration is therefore limited to a simple 2x2 mixing problem: 

(?:)(;)-*(:)■ 

where e% = E(\v)) — E(I max ) and ei = E(\ir)) — E(I rnax ). The problem can be 
solved analytically in two different ways. The first method (A) determines V based 
on strict requirement that "spurious" solution should be placed at zero energy 
Ai = 0. This corresponds to V — y/eiei and A2 = e\ + £2- This is energy of physical 
Imax — 1 state relative to energy of I max state. The results of this procedure for the 
SkO-SHF calculations are shown in the right hand side of Fig. |4] (black triangles). 




_» f-"- 1 ^ 1 ±1 1 l: 1 j i ■ x± 1 ±1 1 l: . 

42 Sc 44 Sc ^Ti 46 Ti 47 V 42 Sc 44 Sc 44 Ti 46 Ti 47 V 



Fig. 4. Energy difference between Imax and Imax — 1 states for /"/ 2 configuration. Left panel 
shows the SkO-SHF solutions |7r) and \u) which break rotational symmetry. Right panel show 
the effect of symmetry restoration according to the methods A and B, respectively. In both cases 
results are compared with experimental data (dots) and the SM calculations (circles). 



The method B uses values of a and b coefficients calculated by acting on 

lax 7 Imax 

we obtain: 



\Imax, Imax) state with /_ operator, see Ref. ^2 for details. For example, in 43 Sc 



I- \Imax, Imax) = \Imax, Imax ~ 1) = 2>/3\tv) + V7\u). (17) 

Knowing a and b one may find from (|16[) the value of interaction V as well as 
energies of "spurious" and physical solutions Ai and A2: 

V= r -^1, (18) 



1 - 




ei - 




1 - 




e2 - 


2 

r z ei 



Ai - ^, (19) 



A 2 = i£_L - , (20) 
1 — r z 
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where r = b/a. Physical solution A2 is shown the right hand side panel of Fig. [4] 
Note, that in this case the energy, Ai, of "spurious" state is not equal zero, but 
for SkO force it doesn't exceed ±0.1 MeV (with exception of 45 Ti). Note also that 
method B doesn't work for N = Z nuclei where r = 1. It is clearly seen from 
Fig H that our simple symmetry restoring schemes provide very accurate results. 
More detailed discussion concerning this issue including results for Sly4 force can 
be found in Ref. 
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Fig. 5. Energy difference between I m ax and the lowest Imax — 1 terminating states for d^^f™^ 1 
configuration. Triangles and grey dots denote SM results and experimental data respectively. Black 
dots represent bare SHF solutions. Circles represent symmetry restored SHF results according to 
the method A. Left (right) panel illustrates the SLy4-SHF (SkO-SHF) results, respectively. In 
order to facilitate comparison of isotopic and isotonic dependence of the Sly4 results we include 
additional curve showing symmetry restored SLy4-SHF solution shifted arbitrarily by 455 keV. 
The method B, which is not shown for the reasons of clarity, gives larger correction than the 
method A by ~100keV for both Sly4 and SkO forces. 



For rf^/™/^ 1 configuration there are three possibilities of creating I max — 1 
states at the mean-field level. One may change the signature of either proton or 
neutron in / 7 / 2 sub-shell. These states will be labeled as \tt) and respectively. 
Alternatively, one may change the signature of the proton in sub-shell. This 
Slater determinant will be denoted as Iff). Simple M-scheme counting method shows 
that our Hilbert space contains two-fold I m ax — 1 representation. Hence, to restore 
broken symmetry we have to deal here with 3x3 eigen-problem having two physical 
and one " spurious" \I m ax, Imax — 1) solutions. The most general mixing Hamiltonian 
matrix can in this case be written as: 



V12 V£ 












) 




-" 













(21) 



where ej = E[i] — E[I max ], and i = \tt), \ v), |tt). Two different methods of restoring 
symmetry (called A and B) are proposed below. In both methods the energy of 
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" spurious" solution is rigorously set to zero. The method A assumes additionally 
real interaction and precise knowledge of mixing coefficients a, b and c. The values 
of a, b and c are calculated again by using simple angular momentum algebra, 
in particular by applying J_ to \I ma x, Imax) reference state, see EH for details. 
Method B admits complex interaction and sets no further constraints on a, b and 
c coefficients. 

In both cases the problem can be solved analytically. Method A gives: 

-eia 2 - e 2 b 2 + e 3 c 2 



Vv. 



2ab 

-era 2 + e 2 b 2 - e 3 c 2 



Vis = V — , (22) 

lac 



V 23 

while in method B we obtain: 



era 2 — e 2 b 2 — e 3 c 2 
2bc ' 



Vn = V e i e 2, V13 = y/eie 3l V 23 = y/e 2 e 3 , (23) 
if ei, e 2l e 3 are of the same sign or: 



Via = V-eiea, Vi 3 = -i\/-eie 3 , V23 = V e 2 e 3> (24) 

for ei < and e2,e 3 > 0. Other possibilities are analogous. Knowing the values of 
interaction one may find the energies of physical I ma x — 1 states. They are equal: 

A± = i (E ± v/A) (25) 



where 



A = X 2 - 4Z, 



ei + e 2 + e 3 , 



Z = eie 2 + eie 3 + e 2 e 3 - |Vi 2 | 2 - |Vi 3 | 2 - | V23 1 2 - 

The results for the lowest I max — 1 states are shown in Fig. [3J One can see, that 
now the effect of symmetry restoration is relatively small as compared to clearly 
dominant mean-field splitting. There is also an interesting difference between Sly4 
and SkO forces. For the Sly4 force the isotonic and isotopic dependence of E(I max ) — 
E(Imax-l) is very well reproduced but there is a constant offset between theoretical 
results and the data. This offset, after symmetry restoration, equals to ~ 455 keV 
for the method A and to ~ 385 keV for the method B. The opposite is true in the 
SkO case. This force does not reproduce details of isotopic and isotonic dependence, 
but reproduces quite well mean value of E(I max ) — E(I max — 1). 



4. Summary 

We have shown that terminating states in A~44 mass region are excellent play- 
ground for testing and constraining various aspects of the nuclear EDF. In par- 
ticular, unification of otherwise quite random coupling constants connected with 
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spin-fields to experimental data and slight reduction of SO strength allow to reach 
very good agreement with experimental data for maximum-spin I max terminating 
states. Our study reveals that the SHF method cannot be directly applied to unfa- 
vored signature I ma x — 1 states. For these states the SHF solutions manifestly break 
rotational invariance at almost spherical shape. After identification of the underly- 
ing SSB mechanism we propose analytical symmetry restoration schemes for both 
[/Jyj] Imam-i an d [^3/2/7/2 ] J m«»-i configurations. It is shown that for [fy 2 ]i ma x-i 
the configuration mixing is absolutely necessary in order to reproduce empirical 
E{I ma x) - E(I max - 1) splitting. For [d^f^hma^-i configuration on the other 
hand the splitting E(I max ) — E{I max — 1) is dominated by mean-field itself and the 
symmetry restoration effect is relatively less important. 
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